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Summary 

Background: Effective countermeasures against emerging infectious diseases require an 
understanding of transmission rate and basic reproduction number (Ro), The Ro for severe acute 
respiratory syndrome (SARS) is generally considered to be >1, whereas that for Middle East 
respiratory syndrome (MERS) is considered to be <1. However, this does not explain the large- 
scale outbreaks of MERS that occurred in Kingdom of Saudi Arabia (KSA) and South Korean 
hospitals. 

Aim: To estimate Ry in nosocomial outbreaks of MERS. 

Methods: Ro was estimated using the incidence decay with an exponential adjustment model. 
The KSA and Korean outbreaks were compared using a line listing of MERS cases compiled using 
publicly available sources. Serial intervals to estimate Rp were assumed to be 6-8 days. Study 
parameters (Rp and countermeasures [d]) were estimated by fitting a model to the cumulative 
incidence epidemic curves using Matlab. 

Findings: The estimated Rp in Korea was 3.9 in the best-fit model, with a serial interval of 6 days. 
The first outbreak cluster in a Pyeongtaek hospital had an Ro of 4.04, and the largest outbreak 
cluster in a Samsung hospital had an Ro of 5.0. Assuming a 6-day serial interval, the KSA 
outbreaks in Jeddah and Riyadh had Ro values of 3.9 and 1.9, respectively. 

Conclusion: The Ro for the nosocomial MERS outbreaks in KSA and South Korea was estimated 
to be in the range of 2-5, which is significantly higher than the previous estimate of <1. 


Therefore, more comprehensive countermeasures are needed to address these infections. 
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Introduction 

The emergence of infectious diseases associated with Middle East respiratory syndrome (MERS), 
severe acute respiratory syndrome (SARS), and Ebola has created unprecedented public health 
challenges. These challenges are complicated by the lack of basic epidemiological data, which 
makes it difficult to predict epidemics. Thus, it is important to quantify actual outbreaks as 
novel infectious diseases emerge. Disease severity and rate of transmission can be predicted by 
mathematical models using the basic reproduction number (Ro)." For example, Ro has been 
extensively used to assess pathogen transmissibility, outbreak severity, and epidemiological 


control.?* 


In previous studies, the Ro for MERS has ranged from 0.42 to 0.92,° which suggests that the 
MERS coronavirus (MERS-CoV) has limited transmissibility. However, these studies typically 
considered community-acquired MERS infections. In this context, nosocomial infections can 
exhibit different reproduction numbers, as the transmission routes for community-acquired and 
nosocomial infections often differ.’ Recent studies have also examined large healthcare- 
associated outbreaks of MERS-CoV infection in Jeddah and Riyadh within the Kingdom of Saudi 
Arabia (KSA). One study reported higher healthcare-acquired Ro values than those from 
community-acquired infections when using the incidence decay with exponential adjustment 


(IDEA) model, which yielded values of 3.5-6.7 in Jeddah and 2.0-2.8 in Riyadh.”” The IDEA 
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model is simple because it does not consider the population-level immune status, which makes 
it especially useful for modelling emerging infectious diseases in resource-limited settings. 

The MERS outbreak in South Korea was associated with hospital-acquired infections. At that 
time, the Korea Centre for Disease Control and Prevention (KCDC) assumed that the outbreak 
had an Ro<1. Thus, the initial countermeasures were not sufficiently aggressive to prevent the 
spread of MERS-CoV infection to other hospitals. Therefore, we used the IDEA model to 
evaluate and compare the MERS Rovalues from the outbreaks in both the KSA and South Korean 


hospitals. 
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Methods 

Data source 

The KSA data were obtained using a line listing of MERS-CoV cases that was maintained by 
Andrew Rambaut (updated on 19 August 2015). The line listing was created using data from the 
KSA Ministry of Health and World Health Organization reports (WHO).”° Since only 44% of the 
cases in the KSA listing included the onset date, hospitalization dates or reported dates were 
used instead. The Korean data were obtained from the KCDC. Among the 186 MERS cases, 178 
had confirmed onset dates. The eight cases with unknown dates of onset were assigned dates 
based on those of laboratory confirmations. All cases in the KSA and Korea were confirmed 
based on laboratory findings. Study parameters (Ro and countermeasures [d]) were estimated 
by fitting a model to the cumulative incidence epidemic curves using Matlab software 


(Mathworks, Natick, MA, USA). 


The data were narrowed down to only the hospital infection cases. Cases with unknown 
transmissions were considered to be hospital infections if a) the patient was in contact with a 
healthcare worker and/or hospitalized patients, or b) the patient was a healthcare worker. Cases 
were excluded if they could not be verified as hospital infections (e.g., zoonotic transmission, 


family contact, or community infection). 
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Model 
We used the IDEA model to estimate the Roas reported previously,’ together with publicly 
available data. The IDEA model is based on the concept that the number of incident cases (J) in 


an epidemic generation (t) that can be counted as: 


I(t) = Ro": (1) 


When an outbreak occurs, epidemic control measures can be implemented, which can, in turn, 
change the Ro. Therefore, the relationship between / and Ro with countermeasures (d) is defined 


as follows: 


I(t) = eo (2) 
(1+d)é 

The Ro and d parameters are estimated by fitting J from model (2) to the observed cumulative 
incidence data of MERS using the least-squares data-fitting method. Since the IDEA model is 
parameterized using epidemic generation time, in the present study, incidence case counts were 
aggregated at serial intervals of 6, 7, and 8 days.”° 
We considered two large outbreaks in each country studied: the outbreaks in Riyadh and 
Jeddah for the KSA, and those in Pyeongtaek St. Mary’s Hospital, and Samsung Seoul Hospital 
for South Korea. The term resnorm is defined as the norm of the residual, which is the squared 
2-norm of the residual; it measures the difference between observed data and the fitted value 


provided by a model. However, since residuals can be positive or negative, a sum of residuals is 
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not a good measure of overall error in the fit. Therefore, a better measure of error is the sum of 


the squared residuals (E), which is calculated as follows: 


E = (F(x, xdata;) — ydata,)’. 


(3) 


The given input data (xdata), the observed output data, (ydata), and F(x, xdata) are the 
functions we wanted to fit, where xdata was an epidemic generation, ydata was the observed 
cumulative incidence data, and F(x, xdata) was equation (2). 

Since the generation times and the estimated values differ according to serial interval times, the 
resnorm changes accordingly. Therefore, to compare the resnorm with the serial interval time, 


the relative resnorm was defined as follows: 


(F (x xdata;)—ydata;)* 


E= i 


ydata; 


The IDEA model was fitted to the cumulative South Korean MERS-CoV case data from the onset 
date of the first case to the onset date of the last case. The outbreak start date was defined as 
11 May 2015 because that was the symptom onset date for Patient Zero, who was the index 
case and caused the outbreak in the Pyeongtaek hospital. MERS patient no. 14 caused the 
outbreak at the Samsung hospital, and his symptom onset date was 21 May 2015. The last case 
of the MERS outbreak in South Korea was observed on 4 July 2015. The KSA MERS outbreak 
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model was fitted using the cumulative incidence data from 28 March 2014 to 2 June 2014 in 


Jeddah and from 20 March 2014 to 29 May 2014 in Riyadh. 


Ethical Considerations 
All data used in these analyses were de-identified publicly available data obtained from the 
WHO, the KSA Ministry of Health website, or KCDC datasets. As such, these data were deemed 


to be exempt from institutional review board assessment. 
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Results 

The KSA outbreaks were relatively large, with 180 cases (over the course of 67 days) in Jeddah 
and 142 cases (over the course of 71 days) in Riyadh. The Korean outbreaks involved 186 cases 
(over the course of 55 days), including 36 cases (over the course of 23 days) in the Pyeongtaek 
hospital, and 91 cases (over the course of 45 days) in the Samsung hospital. Most Korean cases 
(180) were hospital acquired, with the exception of four cases acquired by household 
transmission and two cases with unknown modes of transmission. In the KSA, only two cases 
involved confirmed zoonotic transmission, while a large number of unknown transmissions 
(Jeddah: 99 cases; Riyadh: 69 cases) and hospital exposures (Jeddah: 80 cases; Riyadh: 70 cases) 


were observed (Table |). 


The IDEA model was fitted to the daily KSA and Korea MERS-CoV case data according to the 
onset date. Figure 1 displays the cumulative MERS-CoV case data for the 2014 KSA and the 2015 
South Korea MERS outbreaks. Patient Zero’s symptom-onset date was 11 May 2015; however, 
he was admitted to the Pyeongtaek hospital on 15 May 2015. Therefore, the outbreak was 
assumed to start on 15 May 2015 via a simulation of the Pyeongtaek hospital outbreak. The 
outbreak start date for the Samsung hospital was determined to be 25 May 2015, following the 


same logic (Figure 1). 
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Figure 2 shows the results of the 2014 KSA outbreak. Squares (c), circles (©), and asterisks (*) 
represent data aggregation of the number of cases by serial intervals of 6, 7, and 8 days; the 
curves represent model fits for best-fit parameters. Our estimated Ro values for Jeddah and 
Riyadh were in the range of 3.95-6.68 and 1.92—2.52, respectively, using serial intervals of 6-8 
days. The estimated Ro values for the Korea MERS outbreak were 3.96, 4.91, and 5.95 for serial 
intervals of 6, 7, and 8 days, respectively (Figure 3). Since most cases were related to hospital- 
acquired infections, the Rg for each hospital was also considered. The outbreak in the Samsung 
hospital was larger than that in the Pyeongtaek hospital (the first Korean outbreak). The 
Pyeongtaek hospital exhibited best-fit Rp values of 4.04, 4.23, and 4.39 for serial intervals of 6, 7, 
and 8 days, respectively, while the Samsung hospital exhibited greater Ro values of 5.0, 6.8, and 
8.11 for serial intervals of 6, 7, and 8 days, respectively. Figure 3 shows that the IDEA model 
provided well-fitted curves for the cumulative data regarding South Korean MERS symptom- 


onset dates for all cases. 


Although the IDEA model seemed to be appropriate, the original data never precisely fit the 
model. Therefore, the appropriateness of the model was assessed. Error was evaluated using 
the relative resnorm to find the best-fit parameters. The results indicated that the best-fit Ro 
and serial interval values were 4.9 and 7 days for all cases, 4.39 and 8 days for the Pyeongtaek 


hospital, and 5.0 and 6 days for the Samsung hospital, respectively. Countermeasures (termed 
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182 “d”) increased with each serial interval because the daily effort of countermeasures was 
183 aggregated by serial interval. 
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Discussion 

The clusters of MERS-CoV cases in KSA healthcare facilities occurred from late March to late 
May 2014, while the Korean outbreaks occurred from mid-May to early July in 2015. These 
hospital-based outbreaks exhibited characteristics different from those of community-based 


outbreaks (higher Ro values and case fatality rates) 


The estimated Rais a basic epidemiological variable that is important for selecting appropriate 
countermeasure efforts. However, an emerging infectious disease often has an unknown 
epidemiology, making it difficult to mathematically model. Several methods have been 
proposed to address this issue, including the IDEA model. The Richards model can also estimate 
the Ro using the cumulative daily number of cases and the outbreak turning point (or the peak, 
ae In this context, Hsieh used the Richards model to estimate the Ro values for the Korean 
outbreak as 7.0-19.3. Yet, the Richards model does not consider any countermeasures 


implemented during an outbreak; therefore, it can only be used after an outbreak has peaked. 


The present study used the IDEA model to estimate the Rovalues from the MERS outbreaks in 
the KSA and South Korea. The IDEA model exhibited a good fit: the estimated Ro values for South 
Korea were 3.9-8.0, and the best-fit Rj was 4.9 for a serial interval of 7 days. Conversely, the Ro 
values for Riyadh and Jeddah were 1.9—2.5 and 3.9-6.9, respectively, using serial intervals of 6— 


8 days. Majumder et al.’° used the IDEA model and estimated very similar Ro values of 2.0-2.8 
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for Riyadh and 3.5-6.7 for Jeddah, with serial intervals of 6-8 days. However, the estimated Ro 
values from the present study were much higher than the previously reported values of <1 for 
MERS (the threshold for an epidemic).’° Regardless, the Korean government assumed that the 
outbreak had an Ro value of <1 based on the previous research. The initial criterion for 
quarantine, therefore, was limited to cases of “close contacts,” which were defined as people 
who were within 2 metres of a MERS patient for >1 hour.’® These quarantines—established 
using an incorrectly assumed Ro—resulted in more MERS patients and greater hospital-to- 


hospital transmission.’° 


A serial interval is the interval between successive cases of an infectious disease. This time 
period depends on the temporal relationship between the infectiousness of the disease, the 
clinical onset of the source case, and the incubation period of the receiving case.’ As MERS 
becomes infectious with the onset of clinical symptoms, the MERS latency period equals the 
incubation period. Therefore, the shortest serial interval could be the same as the incubation 
period, and the longest serial interval could be the sum of the incubation period and the 
maximum duration of infectiousness. During the Korean MERS outbreak, several super- 
spreading events occurred because the MERS cases were not immediately isolated upon 
presentation of clinical symptoms.’® Thus, these cases contacted susceptible individuals for up 
to 1 week after the onset of their clinical symptoms. However, most MERS cases with laboratory 
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confirmation were isolated immediately after clinical-symptom onset. In this study, since 
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the incubation period was 2-14 days (median: 6 days), the serial interval was slightly longer 
than the incubation period. The IDEA model with several serial intervals (4-12 days) was used 
and found that intervals of 6-8 days provided the best fit. For the KSA data, even though the 
reported date was used instead of the onset date, the Ro was not affected because aggregated 


data by serial intervals was used in the analysis. 


The IDEA model is limited by the fact that the countermeasures term (d) cannot be compared 
with the d of another model. In this context, an increasing d in accordance with increasing serial 
intervals indicates that the countermeasure efforts are increasing. However, the size of d cannot 
be compared between two or more models of different outbreaks. Nevertheless, the strength of 
the IDEA model is its simplicity because the Rovalue can be estimated using only the cumulative 


number of cases according to the serial interval. 


Conclusions 

The estimated Ro values from the KSA outbreaks (Riyadh and Jeddah) ranged from 1.9 to 6.9, 
whereas the estimated values from the South Korean outbreaks ranged from 3.9 to 8.0. Based 
on these findings, it appears that nosocomial MERS-CoV outbreaks in the KSA and South Korea 
had higher Ro values than the previously assumed values of <1. Although community-acquired 
infections are caused by contact, nosocomial infections are caused by a combination of contact 


and aerosol transmission; therefore, Ro values for hospital infections can be higher than those 


15 


245 


246 


for community-acquired infections. Hence, more comprehensive countermeasures are needed 


to address nosocomial MERS infection and prevent its spread. 
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Tables 


Table |. Characteristics of selected MERS outbreaks in Saudi Arabia and South Korea 


Saudi Arabia South Korea 

Pyeongtaek Samsung 
Jeddah Riyadh Total Hospital Hospital 
Outbreak Onset date 28/3/2014 20/3/2014 11/5/2015 15/5/2015 =. 25/5/2015 
Duration (day) 67 71 55 23 45 
No. of cases 180 142 186 36 91 
Exposure Hospital 80° 70° 180 36 88 
Household 0 3 

Zoonotic 1 1 0 

Unknown 99 69 0 
Status’ Healthcare worker 40 8 39 3 15 
Patient 82 20 36 
Family or visitor 63 13 40 
Unknown 140 134 2 0 0 
Date® Onset date 75 66 178 36 85 
Hospitalized date 85 79 186 36 91 
Reported date 180 142 186 36 91 


: Hospital exposure cases included healthcare workers and individuals who were in contact with a healthcare 
worker or hospitalized patients. 


*The status of cases when they were exposed to MERS. 


>The number of cases with information for onset date, hospitalization date, and reported date of MERS. 
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Figures 


Legends 


Figure 1. Epidemic curves of cumulative cases by selected MERS outbreaks in Saudi Arabia and 


South Korea. 


Figure 2. Best-fit R, by serial intervals of MERS in Jeddah and Riyadh, Saudi Arabia, 2014, using 


the IDEA model. 


Figure 3. Best-fit R, by serial intervals of MERS in South Korea, 2015, using the IDEA model. 
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Figure 1. Epidemic curves of cumulative cases by selected MERS outbreaks in Saudi Arabia and 


South Korea. 
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Figure 2. Best-fit R, by serial intervals of MERS in Jeddah and Riyadh, Saudi Arabia, 2014, using the 


IDEA model. 
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Figure 3. Best-fit R, by serial intervals of MERS in South Korea, 2015, using the IDEA model. 


